Effect of decoherence on resonant Cooper-pair tunneling in a voltage-biased 

single-Cooper-pair transistor 
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We analyze how decoherence appears in the I — V characteristics of a voltage-biased single- 
Cooper-pair transistor. Especially the effect on resonant single or several Cooper-pair tunneling is 
studied. We consider both a symmetric and an asymmetric transistor. As a decoherence source 
we use a small resistive impedance (Re[Z(w)] <C Rq = h/4e 2 ) in series with the transistor, which 
provides both thermal and quantum fluctuations of the voltage. Additional decoherence sources 
are quasiparticle tunneling across the Josephson junctions and quantum /-noise caused by spurious 
charge fluctuators nearby the island. The analysis is based on a real-time diagrammatic technique 
which includes Zeno-like effects in the charge transport, where the tunneling is slowed down due 
to strong decoherence. As compared to the Pauli-master-equation treatment of the problem, the 
present results are more consistent with experiments where many of the predicted sharp resonant 
structures are missing or weakened due to decoherence. 

PACS numbers: 03.65.Yz, 73.23.Hk, 74.50.+r, 74.78.Na, 85.25.Cp 



I. INTRODUCTION 

The quantum mechanical effects originating in coher- 
ent tunneling of Cooper pairs in small Josephson junc- 
tions (JJ) have been investigated actively in recent years 
due to their great potential to be used in nanotcchno- 
logical applications in the future 1 ' 2 . A central obstacle 
has been that the effects are easily decohered by an un- 
controlled coupling between the studied system and its 
nearby environment '. However, due to their high sen- 
sitivity to the environmental fluctuations, the small JJs 
can as well be used as probes of physics in mesoscopic 
low temperature devices 4 ' 5 ' 6,7 . 

In this paper we study theoretically how decoherence 
appears in the I — V characteristics of a voltage-biased 
single-Cooper-pair transistor 1 '* (SCPT). The noise in the 
source-drain line of the SCPT is dominated by quantum 
mechanical fluctuations of voltage due to finite transmis- 
sion line impedance, or more generally due to coupling to 
the electromagnetic environment' 1 (EE). The noise in the 
gate (charge) is due to spurious charge fluctuators 10 ' 11 
(CF) in the materials surrounding the island, as the EE 
noise in this line is usually shielded by a small gate ca- 
pacitance. These decoherence sources are accompanied 
by quasiparticle tunneling across the JJs, where a single 
electron tunnels across a JJ with the cost of breaking a 
Cooper pair. In this work we model the SCPT's trans- 
port properties to the second order in the interaction with 
the relevant unperturbed baths describing EE, CF and 
quasiparticles. This leads to the demand that the resis- 
tance associated with the EE and CF has to be small 
(<C Rq = h/Ae 2 ) and that the tunneling resistances of 
the JJs have to be large (> Rq). 

We study the I—V characteristics in the subgap regime 
where resonant tunneling of Cooper pairs is the essential 
conduction mechanism 12 ' 13 ' 14 . The average current in 
this regime has been previously calculated by using the 



Pauli master equation 9 ' 12,13 ' 15 modeling the populations 
of the SCPT's eigenstates. We call this treatment as 
the coherent tunneling model (CTM). Practically every 
tunneling event of one or more Cooper pairs across the 
JJs in the SCPT can be made resonant by properly ad- 
justing the transport and gate voltages. It follows that 
the CTM-current has quite cumbersome features as the 
subgap region is full of peaks with alternating heights 
and widths. Experimentally, however, the first-order 
(single particle) resonances are well visible 11 but only 
few of the higher-order (several particle) resonances have 
been identified 9 ' 14 ' 16 ' 17 . This kind of "wash-out" of the 
higher-order resonances occurs apparently due to strong 
decoherence caused by the environment as compared to 
the widths of the resonances. The results obtained by 
the CTM are valid only for the strongest resonances and 
therefore only qualitative match between this theory and 
experiments can be obtained. 

We model the system using a density matrix approach 
(DMA) based on a Keldysh type real-time diagrammatic 
technique 1,18 in the Born approximation (leading order). 
As compared to the CTM, the method takes into account 
also nondiagonal contributions in the master equation. 
This enables Zeno-like effects 19 ' 20 in the charge trans- 
port, where strong decoherence drives, or continuously 
"measures" , the SCPT into superpositions of eigenstates 
that are more stable under the decoherence. In Ref. [20] 
this was used to analyze the experimental findings of 
Ref. [21] for a double- island system using the EE as the 
relevant source of decoherence. Also the charge tunneling 
nearby the so-called JQP-cycles has widely been analyzed 
by similar methods 22 ' 23 ' 2 4,25,26 considering quasiparticle 
tunneling as a perturbation. This paper extends these 
results in the sense that the method applies to arbitrary 
resonances and includes all the relevant noise sources in 
the voltage-biased SCPT. We show by numerical simula- 
tions that in typical experimental conditions most of the 
higher-order resonances are indeed lost, unless a care- 
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ful filtering of the relevant noise is done. We also show 
that the ohmic environments (EE and CF) leave differ- 
ent traces to the I — V characteristics so that one can, in 
principle, identify the strength of EE and CF separately 
in a specific experimental setup. In the case of asymmet- 
ric SCPT and strong dissipation, the model produces the 
results of Ref. [27] (incoherent tunneling of Cooper pairs 
across the small JJ). In the opposite limit of weak dissipa- 
tion, the model produces the results of Ref. [28] obtained 
by the CTM. 

The paper is organized as follows. In section II we in- 
troduce the model Hamiltonian for the SCPT and discuss 
its properties. In section III we go through the diagram- 
matic technique used in the modeling, the environmen- 
tal contribution and point out the problems that might 
occur in the DMA and how to circumvent them. The 
Appendix is devoted to specifics of this treatment. In 
section IV we show how average properties can be calcu- 
lated by Laplace-transforming and tracing the equation 
of motion. Section V shows and analyzes numerical solu- 
tions of the problem. The conclusions are given in section 
VI. 



II. THE SCPT HAMILTONIAN 

We model Cooper-pair tunneling across the JJs in a 
nonperturbative way so that the calculation automati- 
cally includes Cooper-pair tunneling processes in arbi- 
trary order. A schematic drawing of the SCPT is shown 
in Fig. 1. The starting point is the Hamiltonian 12 



-ffscPT = 



(Qi -Q2 + Q0) 



2C S 
Eji cos(^i) 



V 



E J2 COs((/? 2 ), 



(1) 



where Qi is the charge tunneled across the i:th J J, a 
conjugated variable to the phase difference ipi, and Eji 
is the related Josephson coupling energy. The capaci- 
tance of the i:th J J is Ci, the capacitance of the gate Co, 
C s = C + Ci + C 2 and Q = C U + (d - C 2 ) V/2 is the 
quasicharge. In the absence of quasiparticle tunneling 
the operator Qi has 2e quantization or, equivalently, the 
cigenfunctions in the ip-space are 27r-periodic 29 . Depend- 
ing whether we are considering an almost symmetric or a 
highly asymmetric SCPT, there are now two convenient 
bases for calculations. 
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FIG. 1: (Color online) A voltage-biased SCPT consisting of 
two JJs (crossed boxes) in series with the voltage source V. 
The island is capacitively coupled to a gate lead. Also is 
drawn the EE described by an impedance Z{uj) and spurious 
charge fluctuators nearby the island (section III). 



Hamiltonian can now be written as 

SCPT ~ ~~ 2Ci 2 Q 



Ej_ 
2 



T^Tq+TqTUtM 



(2) 



where we have defined the charge (Cooper pair) trans- 
lation operators T Q = X) Q \Q + 2e){Q\, T Q = J2q \Q + 
2e)(Q\ and Ej = Ej\ = Ej 2 - One sees that if \j) is 
an eigenstate with an eigenenergy E, then the trans- 
lated state T^' \j) is also an eigenstate with an eigenen- 
ergy E — eV (E + eV). Such states are called equiva- 
lent but belong to different steps of Wannier-Stark lad- 
der called zones. The product state basis \Q)\Q) is the 
most convenient one for a numerical solution in the case 
Ej < Ec = e 2 /2Cs 12 . Since elementary tunneling pro- 
cesses change Qi by 2e, there is no coupling between |0, 0) 
and 0, 2e), where we now use the notation \Q, Q) for the 
product state \Q)\Q). Therefore one can halve the di- 
mension of the Hamiltonian matrix by leaving out the 
state |0,2e), and all the other states obtained from this 
by the elementary tunneling processes. It follows that the 
zones differ from each another by a translation of 4en of 
the feed charge and by —2eVn in the eigenenergy, where 
n is an arbitrary integer. 



B. Asymmetric SCPT and no quasiparticle 
tunneling 



A. Symmetric SCPT and no quasiparticle 
tunneling 

For a symmetric SCPT (identical JJs) the Hamiltonian 
is most easily solved in the island and feed charge basis, 
which is defined as Q = Q\ — Q2 and Q — Q\ + Q 2 - The 



The other special limit is the asymmetric SCPT 27,30 
where by asymmetry we mean that Ej\ Ej2, which 
usually leads to C\ 3> C 2 also. The circuit can be seen as 
a Cooper-pair box 1 (CPB) which is probed, or excited, 
by the smaller J J (probe). A convenient basis for this 
case is defined as Q = Qi — Q 2 , Qt, = Q2, <-P = fx 
</9£ = (fi + <f2- This canonical transformation leads to 



3 



the Hamiltonian 



SCPT 



(Q + Q'o 



-QnV 



2C E 



Eji 
2 



Tn + n 



T^Tq 



(3) 



where C = Co U — CiV (we have assumed that Co <C 
C2). The first two terms on the r. h. s. of Eq. (3) cor- 
respond to the CPB Hamiltonian and the last two terms 
describe single tunneling events of Cooper pairs across 
the probe with simultaneous excitation of the CPB. The 
cigenstates of the CPB are mixed by the probe tunneling 
part mostly when 2eVn, where n is an arbitrary integer, 
matches certain energy level difference of the CPB. 

We solve the eigenstates of the asymmetric SCPT 
representing the Hamiltonian in the product state ba- 
sis |iiQo)|Qs), where \j,Q' ) is the j:th eigenstate of 
the CPB Hamiltonian, which has to be calculated sepa- 
rately, corresponding to the quasichargc Q' . The result- 
ing SCPT eigenstates have the Wannier-Stark structure 
with steps separated by a translation 2en of the charge 
and a change —2eVn in the eigenenergy. 



C. Quasiparticle tunneling 

The preceding treatment is valid when no quasiparti- 
cle tunneling exists, or its effect can be neglected. In the 
opposite case the quantization of the tunneled charge Qi 
must be changed from 2e to e. We treat the quasiparti- 
cle tunneling pcrturbatively causing transitions between 
unperturbed states (section IIIC). Since the unper- 
turbed Hamiltonian i?scPT describes only Cooper-pair- 
tunncling processes, its eigenstates are superpositions of 
states differing 2ne in the tunneled charges, where n is 
an arbitrary integer. In other words the Hamiltonian ma- 
trix is a sum of four blocks each operating in a subspace 
possessing a different parity combination (even or odd) 
of the tunneled charges Qi/ e. The cigenfunctions in each 
of the subspaces can be solved independently by using 
the same Hamiltonian as before, but just properly shift- 
ing the island and/or the feed charge. Only subspaces 
with different parity of the island charge have different 
structures. 



the quantum mechanical noise spectrum is of /-type 10 ' 11 
(and we deal with zero temperature). These decoher- 
ence sources are accompanied by quasiparticle tunneling 
across the JJs. We include the interaction with the rel- 
evant baths by using a real-time Keldysh diagrammatic 
technique 18 ' 32 ' 33 . This leads to a master equation for the 
SCPT's density matrix as the environment is traced out 
of the equations by using the assumption of unperturbed 
reservoirs. We restrict the analysis to the second-order 
calculation in the interaction with each of the baths since 
the higher-order calculation turns out to be actually less 
suitable for a reliable analysis (Appendix). This is justi- 
fied for the EE (CF) for which -R(cf) Rq and for low 
transparency JJs with tunneling resistances Rt > Rq- 



A. Coupling to the electromagnetic environment 

For the case of a dissipative EE, the phenomenologi- 
cal total Hamiltonian describing a voltage-biased SCPT 
in series with the impedance can be presented in the 
form 4 ' 20 



H 



total 



2C\, 



= #SCPT — Cint — 
Cnt 



E 

1 



2d 



Qint 
2Ci nt 

<PR? 



2Li 



(4) 



where the fluctuations caused by the EE are modeled 
by coupling an infinite number of LC-oscillators to the 
phase difference variable <pn across the impedance. The 
EE couples to the SCPT through the term — Qint<z/Ci n t, 
where the charge q passed through the voltage source 
(relative to the equilibrium charge C; n tV) is a conjugate 
variable to cpR. The /-sum part in Eq. (4) can be shown to 
describe correctly the fluctuations across any dissipative 
impedance by properly choosing the parameters L; and 
C/. The extra charging part, which can be written in 
the form (Qi nt — g) 2 /2C nt , can be understood as the 
capacitive energy of the SCPT seen by the impedance. 
The Kirchhoff rules then lead to the identifications (in 
the case of the EE) 



C2C1 + C1Q2 C2 



c s 



c E 



Q + Qt, 



^int l 



(5) 



III. EFFECT OF THE ENVIRONMENT 

The interaction between the SCPT and its environ- 
ment leads to dissipative (open) quantum mechanics ' 1 
and to a net current across the system. Each dissipation 
mechanism shows up in a characteristic way. The noise in 
the voltage across the SCPT is due to a finite transmis- 
sion line impedance (EE) described by Re[Z(w)], which 
we assume to be a constant R at low frequencies. The 
noise in the gate charge is due to the CF, which also can 
be fully characterized by an analogous resistance i?cF if 



and Cnt = (1/Ci + 1/C 2 )~ 1 = C®f (assuming that 
Co <C C\,C<x). This representation of the fluctuations 
corresponds to series LC-oscillators in parallel with the 
SCPT. Alternatively, the EE could also be modeled as 
parallel LC-oscillators in series". 

In the real-time diagrammatic technique the EE is as- 
sumed to be a large reservoir which is not significantly 
affected by the tunneling in the SCPT and is in a good 
approximation described by the last two terms in Eq. (4). 
In this case the environmental part can be traced out and 
described only via its equilibrium properties given by the 
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quantum fluctuation-dissipation theorem 
a(t) = a(-t)* = (8V{t)5V(0)) 
R f°° hw e~ tuJt , 

where SV = q/C- mt , uj c = l/(RC int ) and [3 = l/k b T. 
The calculation now reduces to rules for creating dia- 
grams and the corresponding generalized transition rates 
(Appendix) . The first order diagrams lead to the master 
equation -0 

p(t) = L oP (t) + f E(i - t')p(t')dt' + L N p(t), (7) 

where Lq = (i / h)[- , Hscpt] is the Liouville opera- 
tor, T,(t) the generalized transition rate and Ln = 
(i/H)[-, Qf nt /2C- m t\ its renormalization. The transition 
rate tensor is given by 

S(t) = a'(i)ii„te Lot Li„ t - ia"(t)L int e Lot M inU (8) 

where a(t) = a'(t) + ia"(t), L- mi = (i/h)[-,Q int ] and 
■Mint = Qint] + - The renormalization Ln can be 

shown to cancel certain terms in Eq. (8) (Appendix) and 
to make the renormalized transition rate £(i — t') = £(i— 
t') + 5{t — t')Lff invariant under a 4e translation of Q, or 
2e translation of Qs m the asymmetric case. 

B. Coupling to the charge fluctuators 

The CF can be described similarly as the EE since we 
assume its noise to be of /-type. This means that its 
noise spectral density is proportional to the frequency 
and that Tqp = (temperature of the CF). Also certain 
effective nonvanishing value of Tqf could be used for de- 
scribing dephasing due to thermal fluctuations of the CF. 
However, in this case the model would not be "universal" 
since an ensemble of harmonic oscillators produces differ- 
ent kind of thermal noise as, for example, an ensemble 
of two-level systems. In the description the island charge 
Q = Qi — Q2 is identified as Qi n t of the CF. It is coupled 
linearly to the fluctuating operator q/C(^ , which has the 
desired properties i. e. ohmic dissipation characterized by 
i?CF and a cutoff frequency lo c = l/i?cFCint'- Also the 
renormalization term Q 2 /2C^ has to be introduced in 
order to compensate the change in the effective charging 
potential felt by the SCPT due to the embedded CF 31 ' 34 . 
As a result one ends up with the same total Hamilto- 
nian as in Eq. (4). Therefore the master equation (7) 
applies also to the CF but now with the identifications 
Qint = Qhit = Qi T — Tcf, R = Rcf and C- mt = C iat . 
In the simulations we use the value C^F = C?P. 



C. Quasiparticle tunneling 

The quasiparticle tunneling across the JJs is described 
by using the tunneling-Hamiltonian formalism ' . The 



interaction operator between the quasiparticle reservoirs 
and the SCPT is of the form 

Hn = Tq^ e <g> J2 %M C W + kc -< ( 9 ) 

p£L,k£R 

where is the electron annihilation (creation) operator 
of the state j in the left (L) or right (R) hand side lead 
(quasiparticle reservoir) of the i:th JJ, t\ r the related 

tunneling amplitude and Tq^ £ = \Qi + e )(QA- The 
tunneling rates are calculated using the second-order di- 
agrammatic treatment and dropping out the terms that 
lead to 2e-tunneling of charge (Josephson effect). The en- 
vironmental contribution is now traced in a correlation 
function 

a(*r qp = E I4,»I 2 /(W - f{E k )Y^-^\ (10) 

k,l 

where the sum over the states k and I is done in the equiv- 
alent semiconductor picture 37 and f(E) is the Fermi- 
distribution function. No cut-off has to be introduced 
since we include only the real part of the Laplace trans- 
formed correlation function (Appendix), as the imag- 
inary part would only lead to a small correction on 
the island's capacitance in the case of low transparency 
JJs 32 ' 38 . The factor a(t) l qp can be related to the well 
known superconductor-superconductor quasiparticle cur- 
rent 37 and the operators Tq u£ give the dependence on the 
tunneling direction via the corresponding diagram rules 
(section III CI). 

According to the model above, the equilibrium I — V 
characteristics are always e-periodic in Qq, even with 
zero subgap conductance, and it seems that quasipar- 
ticles can never be neglected. However, the average time 
between the switchings of different parity states of the 
island charge via higher order processes can become very 
long. In Refs. [15,39] the decay of the odd parity state 
(unpaired electron tunnels off the island) was introduced 
to restore the 2e-periodicity in this situation. We in- 
clude the effect by introducing (on the top of the ordinary 
quasiparticle tunneling rate) a constant rate 

7c SC « - * 0(E a - E m ), (11) 
Ze^RrpNi 

for the tunneling of the unpaired electron across the «:th 
JJ 15 . Here Nj is the island's density of states, 0{x) the 
step function, and \a) and \m) are the initial (odd parity) 
and final states of the corresponding Lindblad-cquation 
(section III CI). Theoretically the rate should be en- 
hanced when E a > E m due to singular quasiparticle den- 
sity of states, but we neglect this effect. Throughout the 
paper we use the value 1/Ni = 10~ 2 peV. 

1. The Lindblad form of quasiparticle tunneling 

In the case of quasiparticle tunneling (9) and in some 
parameter region, the numerical solution for the steady 
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state of the reduced density matrix fails by leading to 
negative probabilities for occupancies [this can be stud- 
ied by solving Eq. (22) for k = in the limit s — > 0]. 
In these regions the calculation of the average current 
naturally fails also. This physical inconsistence is possi- 
ble since our truncated (and therefore more or less phe- 
nomenological 40 ) equation is not of Lindblad form. To 
circumvent the problem a proper Lindblad approxima- 
tion of the master equation, which does not change the 
results in the "safe" regions too much, should be used in 
the case of quasiparticle tunneling. 

The Lindblad form can be obtained by two reasonable 
approximations. The first one is the so-called Markovian 
assumption, which in general says that the transition rate 
at time t does not depend on the history p(t'), t' < t. In 
this context one assumes that 



E(t - t')p(t') « 
S(t - t') exp[iH (t - t')]p(t) exp[-iH (t - t') 



(12) 



i. e. the density matrix is approximately a constant in 
the interaction picture on the timescale where the cor- 
relation functions die out. In the diagram language this 
leads to one extra rule (Appendix) which makes the first 
order equations equivalent to the so-called Born-Markov 
(BM) equations 51 in the limit to — > — oo. As the envi- 
ronment is not dependent on the evolution of SCPT, the 
Markov approximation is valid since the memory decay 
time is small compared to the timescale of transitions or 
oscillations of the superpositions 24 . 

In the second step one redefines the nondiagonal tran- 
sition rates in the Lindblad form. The positive-direction 
quasiparticle transition rates across the i:th JJ are in the 
BM-approximation 



- ma \ nb) 



F l (E a — E m ) T l (E b — E n ) 



(13) 



where T 1 ma is the matrix element of the operator TQ ijf , 
and T l (E) = I q - q (E/e)/e is the classical quasiparticle 
transition rate across the i:th J J. The negative direction 
tunneling is obtained by the substitution TQ i E — ► Tq. e . 
The related decay rates (the second-type diagrams) are 



-\Y.^r T L^{E a -E v )8 b 



+T l nv {T l vb yr(E ~E v )8 a 



(14) 



A suitable Lindblad form of the density matrix equa- 
tion is then a sum of four independent operations 



EP = £ 



where is chosen to be 

U it+ =Y, Tl ma\/nEa-E m )\m){a\ 



(15) 
(16) 



and the operator Ui t - is obtained by the substitution 
TQi, e — * Tq e . This gives the same diagonal elements as 
the sum of contributions (13) and (14) but redefines, for 
example, the nondiagonal rates of Eq. (13) to be 



Sa- 
fe- 



E m )^{E b -E n ). 



(17) 



Now, if using E, the transition rates that include forbid- 
den diagonal transitions (i. e. for which E a — E a / < 2A) 
are lost, but other processes are present with approxi- 
mately the same rate due to the high gap energy A. In 
more detail, by assuming that E a — E m > E b — E n > 2A, 
the relevant rates satisfy T(E a — E m )/T(E b — E n ) w 
(E a — E m )/(Eb — E n ) <C 10, from which it follows that 
the Eqs. (13) and (17) are almost equivalent. Numeri- 
cal simulations show that in the parameter region where 
the density matrix calculated from the original master 
equation does not lose its positivity, the above Lindblad 
form leads approximately to the same I — V character- 
istics. It then removes the negativity of the stationary 
state solutions in the problematic regions. 



2. Quasiparticle tunneling thresholds 

There exists a discontinuous jump in the theoretical 
quasiparticle tunneling rates as a function of the volt- 
age since the subgap current occurs only via thermally 
excited quasiparticles. The BCS subgap resistance is ap- 
proximately a constant 12 R q = RTCxp(A/k b T — 1.76), 
which means that the current is negliblc for typical pa- 
rameter values. However, in experiments the thresholds 
broaden to smooth steps and a finite subgap resistance 
persists, probably due to material imperfections 41 ' 42 . 
One extra reason for the threshold broadening in the 
SCPT circuit is the thermal fluctuations of the voltage 
(due to the EE). To include this effect one can proceed 
similarly as in the P(_E)-theory . 

The probability function [P(P)-function] describing in- 
coherent quasiparticle tunneling in the case of a single J J, 
small environmental resistance R, and finite temperature 
T is at low energies (the inhomogeneity of the integral 
equation in Ref. [43]) 



Pe(E) 



tv H + 4£ 2 ' 



(18) 



where L e = A-Kk b TR/R K and R K = h/e 2 . The P e (E)- 
function is damped by the Boltzmann factor e.~ E / kbT for 
larger negative values of E. The EE absorbs (E > 0) 
or emits (E < 0) energy E with this probability in the 
tunneling process. For Cooper-pair tunneling one sim- 
ply has to change Rk to Rq = RkI^, which then leads 
to a width F2 e = 4r e of the P2 e (-E)-function. Gener- 
ally the tunneling of the charge ne leads to the P ne (E)- 
function which is obtained by replacing Rk by Rk(^/ti) 2 



and therefore has the width T r] 



T P n 2 
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To utilize this for a SCPT, we notice that the charge 
seen by the impedance is {C 2 Q[ cal + C 1 Q 1 2 cal )/C s (where 
Qi-cai j g ^ ne charge f th e capacitor i). The change of 
this variable, using the relevant initial and final states, 
determines the amount of charge transferred across the 
impedance in the decay process (in order to regain the 
voltage balance), which then gives the correct value for 
n. For example, in the case of an ordinary quasiparti- 
cle tunneling in a symmetric SCPT, it changes by e/2. 
For the asymmetric SCPT (Eji > E J2 ) but still with 
identical capacitances {C\ = C 2 ) and Eji 3> Eq, the 
incoherent Cooper-pair tunneling across the probe leads 
to a change 2e since the CPB also immediately balances 
the voltage across the larger JJ to zero on the average 
(the energy of the final state has no quasicharge depen- 
dence). In the limit Ej\ <C Eq this does not occur and 
the charge changes by e. The P(-E)-theory gives then for 
the "corrected" quasiparticle tunneling rates [to be used 
in Eq. (15)] 

/OO 
P ne {E')T l {E a -E m -E')dE', 
-OO 

(19) 

where n depends on the states |o) and \b) and on the 
nature of the process. This changes the thresholds from 
sharp to smooth steps. 

IV. CALCULATION OF THE CURRENT 

In studying the average properties, one can simplify 
the density matrix analysis by Laplace-transforming the 
equations (J 4 °° dte~ s ^ t ~ t ° s> ) and considering the limit s — > 

0. The transformation changes the master equation (7) 
to an algebraic equation 

S p(s)-p t0 =£(8)p(s). (20) 

Here S(t) is the (renormalizcd) transition rate tensor, 
which can be expressed as a matrix if p is represented 
as a vector, and p ta is the density matrix at the initial 
time to. For each eigenstate of the SCPT we define index 

1, which gives the number (including sign) of 4e transla- 
tions (or 2e in the asymmetric case) the state differs from 
its equivalent state in the central zone. This defines an 
operator L and the net current can be expressed as the 
time average of 2e-4(Z/(i)}. After the Laplace transform 
this becomes 

I = 2elim s 2 (L(s)). (21) 

Since the transition rates are translationally invariant 
in I (Appendix), we trace the master equation with re- 
spect to this variable. Labeling the states as (a, I), we 
include only diagonal entries in I to the master equation 
(i. e. the elements P(a.h)(b,i 2 ) arc n °t taken into account 
if li ^ l 2 ) to be able to do the trace. Since a state inside 
the central zone can be freely selected to be any of the 



equivalent states, the physics lost in this approximation 
is usually minimized by minimizing the energy differences 
of the states inside a zone. This selection is done sepa- 
rately for each V. The following method can be extended 
to include also nondiagonal entries of the density matrix 
between neighboring zones, etc., but according to our nu- 
merical simulations the results arc not usually changed 
if the zones are chosen as mentioned. An exception is 
mentioned in section V. 

Using the translational invariance we can now do the 
trace over the zones in a Fourier-transform fashion by 
summing the different zone contributions as 

SPmn(k, s) — Pm n (k) 

= I]Nm,*/)(n,/,)( s ) " P^lfKn^)^ 1 '' 
If 

a,b,li ,1 f 

= t^ n (k, S )pa, b (k,s), (22) 
where we have defined the traced transition rate operator 

^(k,s)^^:V(Z^^ e ~ lkl - ( 23 ) 
i 

By inverting the matrix s — S(fc, s), the net current can 
be expressed via the relation 

(L(s))=iTr{±p(k,s)}\ k = 
= ^Tr{[ S -S(fc, S )rV°(M)} |*=o, (24) 
and Eq. (21). 

V. RESULTS 

Qualitatively, depending on the parameters of the sys- 
tem, the I — V characteristics can be divided into three 
regions depending on what kind of processes are domi- 
nant. At low voltages eV < A the characteristics are 
2e-periodic as there are no strong quasiparticle tunnel- 
ing processes. Although higher-order Cooper-pair tun- 
neling can trigger quasiparticle tunneling by releasing 
the needed energy 2A for creating quasiparticle excita- 
tions on both sides of a JJ, this usually occurs at a much 
lower rate than the decay of the odd quasiparticle state 
(section III C). The I — V characteristics are mainly de- 
termined by the coupling of the SCPT to the EE and 
CF. At high voltages eV > 4A (or eV > 2A in the 
highly asymmetric case) ordinary quasiparticle tunnel- 
ing can occur across both the JJs (across the probe) and 
the net current increases steeply almost to the normal 
state value, which is several orders of magnitude higher 
than a typical subgap current. In between is the re- 
gion where, for example, nonrcsonant Cooper-pair tun- 
neling can trigger simultaneous quasiparticle tunneling 
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leading to steps in the I — V characteristics. For applica- 
tions the most important processes in this region are the 
Josephson-quasiparticle (JQP) cycles 8,24 ' 44 , where reso- 
nant Cooper-pair tunneling is accompanied by tunneling 
of quasiparticles. We focus our analysis mainly on the 
higher-order resonances in the region V < A/e. This is 
because the characteristics nearby the usual JQP-cycles 
are not much different from the results obtained previ- 
ously by similar density-matrix approaches 22 ' 23 ' 24 ' 25 ' 26 , 
as the strong quasiparticle tunneling is usually behind 
the broadening of the JQP-cycles overshadowing the ef- 
fect of the EE and CF (T CF = 0). 



A. Asymmetric SCPT interacting with EE 

We first review the results obtained by the CTM for 
the asymmetric SCPT 2, ' 2l \ When represented in the en- 
ergy eigenbasis of the SCPT. the CTM includes transi- 
tions between the diagonal elements of the density ma- 
trix (i. e. populations), but neglects the off-diagonal 
ones. The suitable form of the SCPT Hamiltonian in 
the asymmetric situation, Eq. (3), consists of two parts; 
the CPB Hamiltonian and the probe tunneling operator. 
The name probe is justified since the current I(V) has 
peaks (at least) whenever the energy difference between 
the ground and an excited state of the CPB is equal to 
2eVn, the energy released in the tunneling of n Cooper- 
pairs across the probe. At such a resonance the probe 
tunneling operator mixes the degenerate states |0, Qd)|0) 
and \j, Q' )\2en) (the basis is defined in section II B) and 
causes avoided crossings such that the minimum level 
separation of the corresponding SCPT eigenstates is 

SE = E J2 \(j,Q' \Tl\0,Q' )\, (25) 

for the first order resonances (n = 1). The state 
\j, Q )|2en) usually decays to the state 0, Q' ) \2en) via 
a photon emission to the EE, which is then degenerate 
with |j, Q )|4en) and so on. Thus the resonance leads to 
efficient charge transport. The decay of the CPB states 
via quasiparticle tunneling across the probe (or the larger 
JJ) is similar but costs an energy 2A so that in resonant 
situations it occurs approximately above V = 2A/3e 
(V = A/e), see section VC. The full width at half max- 
imum of the corresponding I(V) peak is determined by 
the width of the mixing region which for n = 1 is approx- 
imately SE/e. For higher order resonances (n > 1) also 
virtual intermediate states contribute and the splitting 
[and the I(V) peak width] has to be calculated sepa- 
rately. For a more detailed discussion of the CTM in the 
asymmetric situation see Ref. 28. 

Fig. 2 shows the /(V) resonances for an asymmetric 
SCPT with Ejx/Ec = 5 calculated numerically using 
the CTM. The band structure (quasicharge dependence) 
of the eigenstates is evident even though Eji/E c > 1. 
The strong resonance oscillating in between the voltages 
40 fxV and 60 fxV occurs since the energy released in 
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Qo/2e 

FIG. 2: (Color online) A two-dimensional map of the current 
I as a function of Qo and V according to the CTM. A contour 
line for the current 10 pA is drawn for clarity. The resonances 
originate in the first or higher-order Cooper-pair tunneling 
across the probe and simultaneous excitation of the CPB, 
described by I/scpt (3). The resonance positions oscillate 
as a function of Qo with the period of 2e due to the band 
structure of the CPB eigenstates. Since Q' — Qo — C2V (the 
quasicharge in the CPB Hamiltonian), the change in V also 
leads to a change in Q' and therefore the band structure is 
skewed when plotted as a function of Qo- The parameters are 
Ej! = 5E C = 12.5£ J2 = A/2.2 = 100 /ieV, C 2 /d = 0.21, 
T — 50 mK and R = 50 Q. 



a tunneling of a single Cooper pair across the probe 
matches the CPB excitation energy to the second (first 
excited) band. A similar resonance between the volt- 
ages 70 fiV and 120 [iV occurs due to the CPB exci- 
tation to the third band. Their Qo dependence is in 
the opposite "phase" because of the band structure. A 
second-order transition to the third band is located at 
40 fiV for Qo = 0. In the region V G 50 - 70 fiV and 
Qo/2e G 0.4 — 0.8 also first-order transitions from the 
second to higher bands are distinguishable since the sec- 
ond band is populated due to the simultaneous 1 — > 2 
resonance. A second-order transition to the second band 
is at 28 [iV for Qo = and a third-order transition to 
the third band at 26 fxV for Qo = 0. The quasiparti- 
cle tunneling at these voltages is very small and can be 
neglected. 

Fig. 3 shows the I(V) resonances for the same sys- 
tem but now calculated using the DMA. The character- 
istics are similar, but differ mostly for the weak higher- 
order resonances. Especially the resonances between the 
first and the third band have vanished or arc signifi- 
cantly weakened in the range from 20 fiV to 50 fiV. 
The second-order resonance to the first band still per- 
sists but its magnitude is reduced. The peaks due to the 
first-order resonances are of equal height in both mod- 
els. According to our analysis the wash-out of the higher 
order resonances occurs due to two decohering mecha- 
nisms that can be treated independently. The first one is 
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Qo/2e 

FIG. 3: (Color online) A two-dimensional map of the current 
calculated using the DMA and the same parameters as in 
Fig. 2. A contour line for the current 10 pA is drawn. As 
compared to the results obtained by the CTM (Fig. 2), the 
characteristics are similar but differ for the weak higher-order 
resonances. Especially, the higher-order resonances between 
the first and the third band are barely visible. The first- 
order resonances are approximately of equal strength in both 
models. 



the level broadening due to a finite lifetime of the CPB 
excited bands caused by quantum noise of the EE. The 
second one is thermal fluctuations of the voltage across 
the SCPT (thermal noise). 

Fig. 4 shows quantitatively the effect of the quantum 
noise by increasing the resistance i?, and at the same time 
the transition rates since they are oc R [see Eq. (Al)], at 
zero temperature when Ej\/Ec ~ 50 (Ref. [27]). As R 
is increased, the transport is first enhanced linearly in 
agreement with the CTM results (not plotted). In this 
region the relaxation is only a small perturbation as com- 
pared to the SCPT eigcnstate splittings and the reduced 
density matrix is almost diagonal in this basis. By further 
increasing R, the widening (due to the finite lifetimes) of 
the CPB excited states starts to overcome these splittings 
and the transport saturates and finally slows down, ap- 
proaching the incoherent limit / oc l/R. This is because 
during the time evolution the fast relaxation damps the 
system to the CPB ground state (which is a superposi- 
tion of the SCPT eigenstates), and the density matrix 
becomes almost diagonal in the CPB basis. This is a 
manifestation of the Zeno or watchdog effect 19,20 . The 
CPB can still evolve to the excited state but now with 
a much lower (incoherent) rate ex l/R. Therefore the 
Cooper-pair tunneling slows down and the net current 
decreases. A similar effect occurs when the splitting of 
the states behind the resonance is lowered by decreasing 
Ej2 at a constant R. 

Fig. 5 shows the effect of thermal noise on the same 
system but now with a fixed R. It demonstrates that 
strong thermal fluctuations also lead to analogous wash- 
out of certain resonances. In Ref. [27] it was derived that 




FIG. 4: (Color online) The effect of the quantum noise ac- 
cording to the DMA. A first-order resonance between the 
ground and the first excited state (which is a level since the 
band width is very small) is located at V ~ 99 fiV and a 
second-order resonance between the ground and the second 
excited state at V ~ 94.5 fiV. With increasing R, the res- 
onant transport first increases oc 7? (coherent regime) and 
finally starts to decrease oc l/R (incoherent regime). The 
width of the resonance in the incoherent regime is determined 
by the lifetime of the CPB excited state" ' and in the coherent 
regime by the splitting of the SCPT eigenstates. The second- 
order transition has much smaller splitting and therefore it 
is more easily washed out by the quantum noise of the EE. 
The parameters are Eji = 48AEc = 500 fieV, E J2 = 3 /xeV, 
C2/C1 =0.1 and T = 0. No quasiparticles are included as 
their effect is small for A > 200 /xeV. 



the (first-order) incoherent Cooper-pair tunneling prob- 
ability function broadens by = Ai^k^TR/ Rq due to 
thermal noise (see also section IIIC2). The DMA pro- 
duces the same result in the limit of incoherent tunneling 
and similar effect occurs also for the higher-order reso- 
nances. However, if the Cooper-pair tunneling in the 
resonance is coherent, the I(V) peak broadens but does 
not significantly weaken. One sees from Eq. (5) that the 
CPB variable Q is protected from the EE's noise by the 
small C2/C5; ratio. This reduction factor does not ap- 
pear in the Qy. variable, which causes dephasing (energy 
level fluctuations) of the tunneled charge states 45 , and 
delivers thermal noise to the system causing broadening. 
In order to see the higher-order effects experimentally, 
besides minimising the environmental temperature, the 
(low frequency) resistance of the voltage line should be 
minimised also, for example, by nearby electrical compo- 
nents. Numerical results show that the factor T 2e should 
not be much larger than the splitting of the relevant 
SCPT eigenstatc behind the resonance if the Cooper-pair 
tunneling is incoherent (broadening due to the quantum 
noise larger than the splitting). 
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FIG. 5: (Color online) The effect of thermal noise on the 
first and second-order I(V) resonances of Fig. 4 with a fixed 
i? = 50 SI. The evolution in the first-order resonance is mainly 
coherent (linear regime in Fig. 4) and the effect of thermal 
noise is fairly small (the area of the peak increases) . Qualita- 
tively this is because the bottleneck of the current is the slow 
relaxation of the populated excited state. The result is sim- 
ilar as obtained using the CTM (not plotted). However, the 
tunneling in the second-order resonance is partially incoher- 
ent and the increase of temperature tends to broaden (with a 
constant area) and finally flush the resonance. Qualitatively 
this is because the bottleneck of the current is the excitation 
of the CPB which is strongly perturbed by thermal fluctua- 
tions of V. In the CTM the second-order resonance widens 
but does not substantially lose its height (the area increases) 
as the temperature is increased. 
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FIG. 6: (Color online) The current nearby two higher-order 
resonances of Fig. 3 for Qo = while changing Ej2 for three 
different environments. The solid lines correspond to Ej2 = 5 
fj,eV, the dashed lines to Ej2 = 10 fieV and the dotted lines 
to Ej2 = 15 fieV. The calculation for the CF (and no EE) 
is made using an effective resistance _Rcf = (C2/Ce) 2 50 fl, 
the calculation for the EE (and no CF) with R = 50 Q, T = 
50 mK but also with T = 0, i. e. when no thermal noise 
is present. The zero temperature current for the EE and 
CF are similar, except that the latter has no nonresonant 
"background" current oc Ej 2 R/V. The thermal noise of EE 
with T — 50 mK broadens the resonances and flushes them 
in the case of small Ej2- 



B. Asymmetric SCPT interacting with CF 

The effect of the quantum /-noise caused by the CF 
is directly linked to the results of the preceding subsec- 
tion. The characteristics for this situation can be ob- 
tained by simply changing Q™, Eq. (5), to = Q 
(section IIIB) and readjusting the environmental resis- 
tivity to a proper value. In Ref. [10] the value Rqf = 3 SI 
was found to match the experimental results and we use 
this in the simulations. Since the current is second order 
in the charge operators, and neglecting the effect of Qs, a 
similar effect is caused by an EE with R = (Cs / C2) 2 Rcf ■ 
This corresponds to R w 100 SI if C 2 /C E « 1/6 (Fig. 2), 
and R = 300 SI if C 2 /C E « 1/10 (Fig. 4). Therefore the 
effect of the CF is usually dominant for a small C2/C1: 
ratio. 

Numerical calculations show that near resonances the 
I — V characteristics are approximately the same as in 
the case of an EE with R = R C f{C^/C 2 ) 2 and T = 0, 
see Fig. 6. This is a consequence of the fact that the 
excited states of the CPB decay only via the operator Q 
which also gives the dominant contributions to the matrix 
elements between the SCPT eigenstates in resonant situ- 
ations. However, the nonresonant current cx Ej 2 R/V 5 ' 28 
is not present in the case of CF since it is related to the 
fluctuations of Qs- This is the main difference between 



the I — V curves and therefore the low bias behaviour 
distinguishes between the two different sources of noise. 

If the CF and EE are present simultaneously, the ther- 
mal noise of the EE can greatly contribute to the rates 
due to the CF. Contrary to the case of quasiparticle tun- 
neling (section III C 2) , this effect is already included by 
summing up two different transition rate contributions 
(from the EE and CF), following from the fact that the 
Cooper pair tunneling is treated nonperturbativcly. The 
thermal noise leads to similar broadening of the reso- 
nances as in the case of the EE, but is now characterized 
by the resistivity of the EE through the corresponding 
broadening factor. Therefore, even if the EE would not 
flush the weak resonances due to a small R or C2/C5] 
(slow relaxation due to the EE), the CF can drive the 
tunneling across the probe incoherent, which then ac- 
companied by moderate thermal fluctuations of EE can 
lead to wash-out of the higher-order resonances. Simi- 
larly, thermal fluctuations of the CF and quantum noise 
due to the EE can lead to the same effect. 



C. Asymmetric SCPT with quasiparticle tunneling 

The quasiparticle tunneling becomes significant in res- 
onant situations when V > 2A/3e for Eji/Ec 3> 1. This 
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is since the decay of a resonant CPB state (releasing en- 
ergy ~ E a — E Q = 2eV) via quasiparticle tunneling across 
the probe (~ eV) releases the total energy ~ ZeV > 2A. 
The tunneling across the larger JJ is more intense but re- 
leases no energy, and occurs therefore approximately at 
voltages V > A/e. These relations are only qualitative as 
the true threshold voltages depend on the details of the 
energy level structure (band structure if Eji/Ec ~ 1) 
and usually are between A/2e < V < 2A/3e. In analogy 
to the EE or CF, if the decay rate due to quasiparticle 
tunneling becomes larger than the splitting of the SCPT 
eigenstates behind the resonance, the Cooper-pair tun- 
neling turns incoherent and charge transport slows down 
(Fig. 4). For voltages well below A/e a more destruc- 
tive effect usually is the simultaneous e-switching of the 
quasicharge Q' , which then leads to an off- resonance sit- 
uation. The system returns to the resonance, for example 
by a decay of the odd parity state of the island, but at 
a much lower rate than typically at resonance. How- 
ever, the switching has no effect in the case of a double- 
resonance point, where the resonance occurs simultane- 
ously in both parity states (for example when Q' = e/2), 
or when the relevant states of the CPB have practically 
no quasicharge dependence {Ej\/Ec ^ 1)- 

Fig. 7 shows the I — V characteristics of the same sys- 
tem as considered in Figs. 2 and 3 but now reducing 
the superconducting gap to A = 150 /LteV and using the 
DMA. At low voltages the characteristics are 2e-periodic 
in Q (as before) but approximately above A/2e, start 
to show e-periodicity as the quasiparticle tunneling in- 
creases. The I— V characteristics due to the (first-order) 
excitation to the third band are now modified in compar- 
ison with Fig. 3. The peak due to the excitation to the 
bottom of the third band is still present but otherwise the 
transition is visible only near the double resonance points 
(pointed by the arrows). The characteristics in this re- 
gion are unchanged when using the CTM (not plotted) 
which indicates that the Cooper-pair tunneling is mainly 
coherent. For A = 100 /ieV the tunneling in this region 
becomes partly incoherent as the quasiparticle tunneling 
across the larger JJ also starts to relax the CPB. In this 
case only the double resonance point is highlighted. 



D. Symmetric SCPT 

Let us now consider resonant tunneling in the sym- 
metric SCPT. We focus on the case Ej < Ec since 
our choice of the basis (section II A) is not optimal for 
numerical simulations in the opposite situation. Also 
the results of the preceding subsections are less valid 
for Ej ^> Ec, where the nonresonant background cur- 
rent (which is cx Ej in the symmetric case) takes a 
more dominant role in the I — V characteristics. Since 
Ej < Ec, the regions of resonant single-Cooper-pair pro- 
cesses can be found analytically by demanding the degen- 
eracy of the charging energies [the first two terms on the 
r. h. s. of Eq. (2)] before and after elementary tunnel- 
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FIG. 7: (Color online) A two-dimensional map of the cur- 
rent calculated using the DMA and the same parameters as 
in Fig. 2, except A = 150 fieV. A contour line for the current 
10 pA is drawn. The resonances below V ~ 80 /iV maintain 
2e-periodicity as the quasiparticle tunneling is very small at 
these voltages. The resonances above this start to show e- 
periodicity as the decay process releases enough energy for a 
Cooper pair to break and the quasiparticle to tunnel across 
the probe. The arrows point to the double-resonance points of 
the third band (see text). All quasiparticle tunneling thresh- 
olds have been broadened by the factor F2 e (section IIIC2). 



ing processes 41 . The process where the charge le tunnels 
across the left JJ and re across the right JJ (both to the 
positive direction) is resonant when 
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where en is the initial island charge. This releases en- 
ergy for larger values of the voltage and if it includes 
quasiparticle tunneling (i. c. I or r is an odd number 
and the process is not the decay of the odd parity state) 
the term 4A/e(Z + r) has to be added to the r. h. s. of 
Eq. (26). After such a quasiparticle process the charge on 
the island changes parity, which effectively means shift- 
ing Qq by e. Efficient charge transport including quasi- 
particles requires cycles of rapid processes in both parity 
states. In the following we label the tunneling processes 
as (I, r, n). The locations of some resonances and thresh- 
olds are drawn in Fig. 8. 

Fig. 9 shows the I—V characteristics for Ej/Ec = 0.3 
and typical EE and CF calculated using the DMA. At low 
voltages and Qq ~ e the current is enhanced as the pro- 
cesses (0, 2, 0) and (2, 0, —2) (see also Fig. 8) with match- 
ing initial and final states produce an efficient charge 
transport cycle. When quasiparticle tunneling is not in- 
cluded, the area of the enhanced current is V-shaped, 
fenced by the resonance lines, and extends above the 
plotted region with slowly decreasing current. In this 
case the shape of the I — V curve (for a fixed Qq) is more 
a step than a Lorcntzian. This reflects that, in contrast 
to the asymmetric case, the single-Cooper-pair tunneling 



11 



eV/E c 



(0,2,0) (2,0,-2) 




FIG. 8: (Color online) Some Cooper-pair tunneling 
resonances (solid lines), quasiparticle-tunneling thresholds 
(dashed lines) and the odd-parity-state decay thresholds (dot- 
ted lines) in a symmetric SCPT with Ej < Ec- The quasi- 
particle thresholds are calculated assuming A = 2E C . 
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FIG. 9: (Color online) A two-dimensional map of the current 
across the symmetric SCPT according to the DMA. Some 
of the resonances and thresholds are identified in Fig. 8. In 
the white regions the current exceeds the plotting range. All 
quasiparticle tunneling thresholds have been broadened by 
the factor F2 e . The parameters are Ec = 3.3Ej = A/2 = 
100 MeV, R = 50 Q, T = 50 mK and R C f = 3 fi. 



and simultaneous photon emission to the EE or CF is a 
strong process above its threshold. Including the quasi- 
particle tunneling changes the characteristics from the 
smooth decrease to steps where quasiparticle tunneling, 
in most cases accompanied by tunneling of one or more 
Cooper pairs, becomes possible. In the step the current 
can also be reduced if the new distribution does not sup- 
port the charge transport processes below the threshold. 

Strong higher-order resonances at low voltages are the 
processes (2, 4, 0) and (4, 2, —2). They are visible in Fig. 9 
as straight lines and are not washed out due to thermal 
fluctuations of the EE since their minimum splitting 1 
« 8lE 3 j/l280E c « 0.7 ^eV is not dominated by the ther- 



mal broadening factor Y2e ~ 0.4 /xeV. Also the quantum 
noise and the quasiparticle tunneling are not too intense, 
except for eV > 3Ec/2 when ordinary quasiparticle tun- 
neling becomes possible for one of the resonant states. 
In Fig. 8 this corresponds to, for example, the cross- 
ing point of the processes (4,2,0) and (0,1,2). Above 
this the CTM current increases but the DMA current 
weakens significantly due to a strong Zcno effect. The 
resonance is also seen in the experiments 9 ' 17 , but disap- 
pearing even before the onset of quasiparticle tunneling, 
perhaps indicating a more pronounced subgap quasipar- 
ticle current. The double-resonance points eV = 2Ec 
and Qo/2e = 0.25 or 0.75 are the locations of the double 
JQP-cycles but for this choice of parameters the ordi- 
nary quasiparticle tunneling is not possible (the points 
are below 2 A — Ec) and the transport occurs via higher- 
order processes. The thresholds above these resonances 
are due to quasiparticle tunneling accompanied by tun- 
neling of two Cooper-pairs. Several low-order processes 
are behind the strongly enhanced current just above the 
double-resonance points. 

Above the range of Fig. 9 at Qo = and eV = 4E C , 
the first-order processes (2, 0, 0) and (0, 2, 0) are resonant 
simultaneously but with nonmatching states. The CTM 
produces a strong current peak through higher-order pro- 
cesses 12,13 but in the DMA the current is only slightly en- 
hanced as the strong quasiparticle current destroys the 
coherence of the weak higher-order process. Since many 
resonances occur simultaneously at this point the DMA 
needs to be done by considering also the "neighbouring" 
nondiagonal states in the zone parameter (section IV) 
and by a careful choice of the states that belong to the 
same zone. The current nearby the ordinary JQP-cycles 
is also reduced due to intense quasiparticle tunneling. 
Actually, if strictly using the Ambcgaokar-BaratofP'' 46 
values for the Josephson coupling energies (with a typi- 
cal energy gap for an aluminum film), this is always the 
case for a symmetric SCPT since the quasiparticle re- 
sistance and the Josephson coupling energy cannot be 
varied independently (for a constant energy gap). This 
limitation does not apply to an asymmetric SCPT . 

Finally, we have studied the effect of the CF alone and 
compared to the case of the EE. Since = Q/2, it 
seems that the EE does not induce quasicharge fluctu- 
ations (tx Q), which mainly cause the relaxation in the 
asymmetric case. Thus qualitative differences could be 
expected in the case of the CF. This is, however, mis- 
leading since this kind of immunity is true only for fully 
symmetric resonances. In the asymmetric basis, which 
can be used in most resonant situations for the evalua- 
tion of the current as well, they exist with a capacitivc 
shielding (C2/C5:) 2 = 1/4. Numerical results show again 
that the CF lead to similar I — V characteristics as the 
EE with R = (C s /C 2 ) 2 Rcf = 4i? C F and T = 0. This re- 
lation tells however that an EE with i? = 50 £1 is qualita- 
tively about four times more destructive than i?cF = 3 O, 
i. e. the noise of the EE usually dominates the noise of 
the CF in the case of symmetric SCPT. Special situations 
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are the regions where similar tunneling processes across 
both of the junctions arc resonant simultaneously, where 
the effect of the CF vanishes. 



VI. CONCLUSION 

We have analyzed typical decoherence mechanisms 
present in mcsoscopic superconducting devices through 
their effect on the I — V characteristics of a voltage- 
biased SCPT. We have shown that each of the envi- 
ronments leave different traces to the characteristics on 
which ground they can be identified. We have also shown 
that the higher-order resonances obtained by the CTM 
tend to be washed out by quantum or thermal noise sup- 
plied by the environment. This explains why only few 
of them have been detected in the experiments. In order 
to see other resonances, the relevant noise sources caus- 
ing the wash-out should be filtered as well as possible. 
Theoretically one can calculate the minimum splitting of 
the SCPT eigenstates causing the resonance and obtain 
analytic relations for preventing the wash-out. In simple 
terms the splitting should not be essentially smaller than 
the decay rate due to the EE, CF or the quasiparticle 
tunneling, or the resonance is not seen due to the Zeno- 
effect. If it is smaller, but not essentially smaller, the 
thermal noise of the EE, described by a similar broaden- 
ing factor, tends to flush the resonance. The effect of the 
CF's thermal noise, which is usually of l//-type, was not 
studied but could also lead to a similar effect. Finally, 
if quasiparticle tunneling exists (but is not too intense), 
the resonance lines are usually only seen near double- 
resonance points, because in other regions the quasipar- 
ticle tunneling switches the system to an off-resonance 
situation slowing down the charge transport. 
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APPENDIX A: CALCULATION OF THE 
TRANSITION RATES 

Here we give rules for calculating the generalized tran- 
sition rates between the SCPT's density matrix entries 
due to the coupling with the EE or CF. The derivation 
of the diagram rules is based on a series expansion of the 
system's time evolution operator in the interaction pic- 
ture 18,32 '' !,! . We also analyze the renormalization, trans- 
lation invariance and validity of the higher-order calcu- 
lation. 



Type 1 Type 2 

t' t f t 




FIG. 10: (Color online) The second-order irreducible di- 
agrams describing interaction with the environment. The 
dashed lines represent the mirror diagrams. Their contri- 
bution is the complex conjugate of the original ones with 
switched a <-> b and m «-> n (a mirror rule). The arrows 
point to the directions of time in the two branches. 



1. Interaction with the EE and CF 

The interaction between the SCPT and the EE/CF is 
described by "environmental points" located at arbitrary 
branches and times in the Kcldysh diagram and are pair- 
wise linked to each other by "environmental lines" . The 
effect of the renormalization term is described by "RN 
pairs" . A RN pair is located at any branch and time but 
the points forming the pair have the same branch and 
time. The contributions from the EE/CF and from the 
renormalization have to be calculated to the same order 
in the sense that one RN pair is equivalent to one environ- 
mental line. The generalized transition rate £^™(i — t') 
is a sum of all irreducible diagrams starting from t' and 
ending to t, meaning that a vertical line between times 
t' and t always cuts an environmental line. The rules 
for constructing the transition rate corresponding to a 
diagram are the following. 

1. Each point or RN pair in the upper branch pro- 
duces a factor — i/H and in the lower branch a factor 

i/h. 

2. Each EE/CF line contributes a factor a(t 2 — t{)^*\ 
where the complex conjugation takes place if the 
left end of the line, corresponding to the smaller 
time t\, is located at the lower branch, and each 
RN pair a factor l/2Ci n t. 

3. Each point produces a factor 

(f\Qba\i)e# i ~ t 'X Uf ~ Ui) where \i) is the enter- 
ing state, \f) the leaving state with respect to 
the direction of the branch, fttOj an eigenenergy of 
the SCPT's eigenstatc \j) and t the timing of the 
point. 

4. An integration over the timings of EE/CF points 
or RN pairs, located between the times t 1 and t, is 
made and the result is multiplied by e^*~* )(<* ; » _a '">). 

The second-order irreducible diagrams describing cou- 
pling with the EE/CF are shown in Fig. 10. From 
these one obtains the transition rates (8). Only one 
RN pair can exist in a second-order graph, with no ac- 
companying EE/CF points. Thus the only irreducible 
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RN diagram is the one which starts and ends at the 
same time t producing the lowest order renormaliza- 
tion operator i/h[-,Qf nt /2C- m t\. Note that according to 
the rules the exponential factor of the first-type solid 
line diagram in Fig. 10, for example, is e^* _t 'f" l "" m ', 
The Markov approximation of the master equation (used 
in section III C 1) would effectively add an extra factor 
e i{t-t ){uj a -uj b ) ^ Q k e transition rates, as one drops out the 
history of the density matrix from the equations. There- 
fore the exponential factor of this diagram in the BM 
approximation is e^ t_t i, e. the energy release 

in the upper branch process (and the lower branch in the 
case of the mirror diagram). We have tested that our 
numerical results are unaffected by the Markov approxi- 
mation in the cases of EE and CF. 



2. Laplace transform and renormalization of the 
transition rates 

The Laplace transform (or time averaging) of the 
second-order transition rates can be reduced to an inte- 
gral of the form / °° dte i{ - u+i8 ^ a(t) , or in the fourth-order 

/ °° dt J*dt 2 J* 2 dhe^+^+^+^ait-t^afo-h), 
from which all the elements of the transformed transition- 
rate tensor can be constructed. At least up to the fourth 
order, the integrals can be evaluated analytically by first 
performing the time integration, then taking the limit 
s — > and, if needed, using the residue theorem. The 
second-order integration gives 

_i^£Re[Z t (a;)] + i^MWM"), (Al) 

Z Z7T 

where 

-2Fe[*(l + i^)], (A2) 

and \&(x) is the digamma function. For this expression we 
have assumed that h/3oj c /2ir = n + 1/2, where n is a large 
integer. The first term on the r. h. s. of Eq. (Al) is the 
usual fluctuation spectrum of the EE or CF and describes 
dccohcrence via emission or absorption of photons. The 
last two terms are purely imaginary. To a good approx- 
imation they produce a constant —iU/2C- m t at usual fre- 
quencies which are well below lo c . They induce coherent 
oscillations and reflect the effective potential caused by 
the EE or CF. In comparison, the renormalization Ln 
in Eq. (7) contributes through the same, but opposite 
factor ih/2C ul t- As this enters E through the same ma- 
trix elements as obtained from the second-type diagrams 



in Fig. 10, there is cancellation of the imaginary terms. 
The constant complex part also vanishes in the sum of 
the first-type diagrams, because the mirror diagrams con- 
tribute by complex conjugated terms. 



3. Validity of the expansion in the case of the EE 

By looking at Eq. (4) one can see a problem: as charge 
starts flowing across the SCPT, starts to increase, 
and as q does not, the charging energy in Eq. (4) starts 
to increase also. So qualitatively speaking, in the exact 
product-state calculation the current should stop after 
some time the decoupling is made, in order to conserve 
the energy. Also the transition rates due to (higher- 
order) processes that include large changes in the feed 
charge should be suspicious since the bath does not "fol- 
low" the state of the system. This is clearly not what 
happens in a real physical situation. The effects follow- 
ing from the assumption of a static EE restricts the usage 
of the real-time diagrammatic technique for this problem. 

The extra imaginary term in Eq. (Al), dependent on 
ty(u>), is not a constant and induces (small) spurious dy- 
namics. This indicates that the effective potential felt by 
the subsystem in this treatment is not exactly Qf nt /2Ci nt . 
However, in the formal solution of Hcisenbcrg equations 
of motion the embedding of the bath (change in the ef- 
fective potential) vanishes exactly . In the second-order 
calculation the term ty(uj) does not contribute in the di- 
agonal terms, as the complex part vanishes, and the con- 
tribution through nondiagonal terms is small. Still dia- 
grams that posses terms like (i\Q^\i), which depend on 
the "position" in the W-S ladder are a source of small 
translational variance, but their effect seems to be small 
as long as |(<3hi?)l/ e ~ 1- The effect can be studied by 
shifting the central zone. 

In the fourth-order calculation the previous change in 
the effective potential and the translational variance seem 
to enhance each other. This occurs because after each in- 
tegration the corresponding RN diagrams do not remove 
the embedding of the EE completely but leave behind the 
discussed extra term. This is then enhanced by terms like 
(i\Qf^\i), which can exist two times in diagrams that de- 
scribe decay, producing strong translational variance and 
dissipation, which is essentially larger than in the first- 
order calculation. Also nondiagonal elements suffer from 
the same effect. If one now removes the corresponding 
spurious terms after each integration from the calculation 
"by hand" , the discussed behaviour vanishes and the final 
correction is small compared to the first order calculation, 
as expected for R/Rq <C 1. So for analysing the effect 
of the EE in the case of larger R, where the higher order 
effects should become dominating, the model cannot be 
used in this form. 
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